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Preliminary results from simulations with 2+1 dynamical quark flavors at a lattice spacing of 0.09 fm are 
combined with earlier results at a = 0.13 fm. We examine the approach to the continuum limit and investigate 
the dependence of the pseudoscalar masses and decay constants as the sea and valence quark masses are separately 
varied. 



The MILC collaboration is engaged in a pro- 
gram of QCD simulations using three flavors of 
dynamical quarks — two light and one strange. 
These simulations use the "af^^" action for stag- 
gered quarks, which removes the tree level order 
discretization errors]^]. Results for scahng of 
hadron masses, for dynamical quark effects on the 
static quark potential, and for the light hadron 
spectrum at a « 0.13 fm have been reported in 
Refs. 1^. We are continuing these simulations at 
smaller quark masses and on finer lattices, with 
a « 0.09 fm. These new runs allow us to ex- 
tend our studies of the approach to the contin- 
uum limit, and to study some new problems such 
as exotic hybrid mesons with staggered quarks. 
Several projects based on lattices from these sim- 
ulations were reported in other talks at this con- 
ference 

Before exploring more technical points, we 
would like to emphasize that the big picture of 
hadronic physics coming from lattice simulations 
is encouraging. For example, in Fig. |] we show 
some hadron masses from our three flavor a « 
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0.13 fm calculations compared to experimental 
results. In this plot the masses (squared masses 
for the pseudoscalars) were simply extrapolated 
linearly to the physical light quark mass, and no 
continuum extrapolation was made. Here the p 
mass was used to set the scale, and the n mass to 
set the light quark mass. (We think that the ai 
and bi energies are small because of couplings to 
two meson states, and we are investigating this 
further.) 

A preliminary indication of the effect of 
nonzero lattice spacing (the continuum extrap- 
olation) is seen in Fig. ||, where the nucleon and 
vector meson (p) masses are plotted in units of 
ri. The leading corrections are expected to be 
proportional to a^g^, or about half the size in the 
0.09 fm runs as in the 0.13 fm runs. One can 
see that the corrections are small, but possibly 
larger in three flavor QCD than in the quenched 
approximation. 

We have used the static quark potential to fix 
the lattice spacing in our simulations. The po- 
tential is also useful in phenomenological models 
of hadrons. In our earlier simulations we showed 
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Figure 1. 2+1 flavor lattice spectrum results 
compared with the real world. These are from 
simulations with a w 0.13 fm (no continuum ex- 
trapolation), with a simple linear extrapolation 
to the physical quark mass. 
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Figure 2. Nucleon and p masses in units of 
ri (upper and lower bands, respectively). The 
open symbols are a = 0.13 fm. Bursts and crosses 
are a = 0.09 fm (preliminary). Fancy plusses are 
physical values, using ri = 0.35 fm. 



clear effects of sea quarks on the shape of this 
potential. With our 0.09 fm runs in progress, we 
can begin to investigate the nonzero lattice spac- 
ing corrections to this shape. Figure || shows one 
shape parameter, rp/ri. (ro and ri are defined 
by r'^F{r) ~ 1.65 and 1.00 respectively.) The 
horizontal axis is basically the dynamical quark 
mass, with the quenched limit at the right and 
the physical point at 0.033. The kink at 0.45 oc- 
curs because for quark masses larger than this we 
varied all three quark masses together, while for 
light quark masses smaller than this the strange 
quark mass was held fixed. The crosses are pre- 
liminary 0.09 fm results. We can see that there is 
a small upward shift as the lattice spacing is low- 
ered, leading to a very crude estimate of a shift 
of 0.02 between a = 0.13 fm and a = 0. 

The relatively small quark masses accessible 
in simulations with staggered quarks allow us 
to clearly see the effects of chiral logarithms. 
As Sharpe and Shoresh have emphasized^, par- 
tially quenched calculations can be used to de- 



termine parameters in the chiral lagrangian that 
describes the interactions of low energy pions. 
Figures § and || show partially quenched pseu- 
doscalar masses and decay constants as a func- 
tion of valence quark mass for several values of 
the sea quark mass. Curvature from chiral loga- 
rithms is clearly visible, as is the effect of varying 
the sea quark mass. Fitting this data will require 
taking into account the effects of the remaining 
flavor symmetry breaking, as discussed by Claude 
Bernard in these proceedings [||. 
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Figure 3. A shape parameter for the static quark 
potential. Octagons and squares are 0.13 fm runs 
for three equal masses and 2+1 flavors respec- 
tively. The crosses are preliminary 0.09 fm re- 
sults. The burst is an extrapolation to the phys- 
ical value, of {rriTr/mpY . 
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Figure 4. Partially quenched Mpg/{mi + 7712) 
for a ~ 0.13 fm. mi and 7712 are valence quark 
masses. Light sea quark masses are 0.2, 0.4 and 
0.6 m,. 



REFERENCES 

1. K. Orginos and D. Toussaint, Phys. Rev. 
D 59:014501 (1999); Nucl. Phys. B (Proc. 
Suppl.) 73, 909 (1999); K. Orginos, R.L. 
Sugar and D. Toussaint, Phys. Rev. D 60 

(1999) 054503, Nucl. Phys. B(Proc. Suppl.) 
83, 878 (2000); G.P. Lepage, Phys. Rev. D 
59, (1999) 074502. 

2. The MILC collaboration: C. Bernard et ai, 
Phys. Rev. D 61, (2000) 111502; Phys. Rev. 
D 62, (2000) 034503; Phys. Rev. D 64 (2001) 
054506. 

3. Work presented by C. Davies, C.E. DeTar, 
Z. Fu, E. Gregory, S. Gottlieb, J. Hetrick, 
S. Tamhankar, and M. Wingate, these pro- 
ceedings 

4. S. Sharpc and N. Shorcsh, Phys. Rev. D 62 

(2000) 094503. 

5. C. Bernard, Phys. Rev. D 65 (2002) 054031; 
also these proceedings. 



0.24 



0.18 



0.16 





' 1 ' ' ' ' 


















□ 








o.em^m^: 


- s° 
" to 




0.4m„m,- 


m,,,= 
0,0, □: mi = 


O.Sm^m^. 
= (tt) - 




1 +,x,*: = 


= (K) : 



0.0 0.1 0.2 0.3 

ri(mi+mp) 
Figure 5. Partially quencned fps- These were 

obtained from point-point propagators, with a ~ 

0.13 fm. The "+" marks the physical /,r and Jk, 

using ri = 0.35 fm. 



